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1. INTRODUCTION 

The recent approach by most utilities in fulfilling the ever-growing electricity demand is distributed 
generation (DG). There are many advantages as a result of DG installation in a distribution system. As has 
been reported in [1-4], the DG is able to improve power quality, voltage stability and reduce power losses 
due to various factors like poor switching operation, voltage dips, interruption, transient, and network 
disturbances at consumer side. It can also improve the efficiency of distribution systems if proper installation 
was carried out. As electric power generation, the DG system can provide real power to the distributed 
system network or the consumer side network [5]. With the purpose of providing a source of real power, 
the DG will not provide any reactive power to the system. However, the reactive power can be provided 
under certain circumstances, especially the one that relates to voltage stability. Due to the power loss which 
will result in insufficient power supply, mostly the location for installing the DG is at consumer side network. 
In addition, the inconsistence of power demand at consumer side is also one of the factors to install a DG. 
The optimal installation, sizing and operation of DG units connected directly to the network can be 
determined by several methods. Mainly, there are four methods can be used in optimal DG placement such as 
the L-index method, Voltage Stability Margin Index (VSMI) method, Voltage Performance Index (VPI) 
method and Loss Sensitivity Factor (LSF) [6-10]. Placement via LSF requires simulation of load flow to 
obtain the resulted real power loss with respect to the DG’s injected real power (i.e. the DG capacity). 
Whichever location in the distribution system that has the lowest LSF value will be noted as the most suitable 
location for DG placement. The power rating for the DG which will be connected to the distribution system 
depends on the capacity of the system. This is divided according to the DG sizes as follows: 1 kW to 5 kW 
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for micro, 5 kW to 5 MW for small, 5 MW to 50 MW for medium and 50 MW to 300 MW for large DG 
size [11]. The technology used in distributed generation can be divided into several categories, such as 
renewable energy generation, modular generation and combined heat and power (CHP) generation [5]. The 
modular generation generates electricity from battery storage and it cannot provide large power demand. 
The combined heat and power (CHP) generation or simply co-generation is the one that generates electricity 
and useful heat from steam using the combination of fossil fuel, engine and gas turbine. Some of these types, 
however, have small disadvantages as they are not renewable and have small effect to the environment. In the 
renewable energy, it generates electricity from natural source such as geothermal, wind, hydro, biomass and 
solar. This type is greener to the world and the energy source is infinite. 

There are many researches on optimal installation of distributed generation using various 
algorithms [12, 13]. Evolutionary Programming (EP), Particle Swarm Optimization (PSO), Genetic 
Algorithm (GA) and Differential Evolution (DE) are the examples. In comparison with traditional 
optimization techniques, meta-heuristic algorithms are more complex as they need many variables and 
parameters to be considered for a good accuracy. Firefly Algorithm (FA), developed by Yang in 2008 [14], is 
a meta-heuristic algorithm inspired from the adapted behavior of firefly swarms to develop optimizing 
function with multiple optima as an algorithm [15-18]. The FA is flexible in terms of implementation and 
problem formulation as compared to other algorithms. In [19], adaptive PSO was applied in the problem of 
DG and capacitor bank allocation. The adaptive version was produced by replacing the inertial weight 
equation in the original PSO’s velocity update equation. Experiment on three different test systems verified 
the effectiveness of the adaptive PSO in solving the DG allocation problem. Later, article [20] proposed 
optimal DG allocation using Whale Optimization Algorithm (WOA). The research considers various types of 
DG resources like photovoltaic cells, fuel cells and kVAR compensators. Through comparison with other 
optimization methods, it was observed that the proposed method yields to the lowest power loss minimization 
and highest voltage profile improvement. A research that employed FA for optimal DG allocation was 
conducted in [21]. The research takes into consideration on voltage profile improvement besides power loss 
minimization in IEEE 33-bus radial network. However, the sizing focused on injected real power by DG 
instead of reactive power. Other researches on optimal DG allocation using various meta-heuristic algorithms 
can be found in [22-25]. 

In this study, the FA acts as an engine to find the optimal size of DG capacity. The aim of the 
optimization is to improve the voltage stability of distribution system, and this will be in terms of voltage 
profiles and real power losses. Prior to conducting the optimization, the identification of suitable locations for 
DG placement will be performed using the LS technique. This means that a rank of buses from the 
distribution system with their corresponding LS indices will be obtained as a result of load flow simulation. 
This paper is organized as follows: Section 2 will describe the proposed methodology and techniques to be 
implemented in this study, Section 3 analyzes the results with some interpretation related to voltage stability, 
and Section 4 concludes the overall findings of this study. 


2. RESEARCH METHOD 

The methodology of this study starts with the philosophy of Firefly Algorithm (FA), followed by the 
proposed problem formulation including the set of variables and constraints to be used in the FA and the 
overall algorithm or process specially designed for this study. 


2.1. Firefly algorithm 

The FA was inspired from fireflies behaviour of using the flash signal to attract the other potential 
fireflies [14, 18]. Based on Figure 1, the behaviour of the meta-heuristic FA is based on the following rule: 
all fireflies are considered unisexual and the attractiveness between them is directly proportional to their 
flashing light intensity or brightness. 





Figure |. Firefly movement 
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Thus, the movement of firefly will be from the one which is less bright to brighter firefly. However, 
as the distance between the two fireflies increases, the brightness of a firefly with respect to the other one 
decreases. If no brighter fireflies nearby, the less bright firefly will move randomly. The FA has two 
important parameters, which are light intensity and the formulation of attractiveness. The equation to 
represent the light intensity is given in (1) 


I(r)=Ive” (1) 


Where, 
I, : light intensity of the source 
I: fixed light absorption 
R : singularity (r = 0) 
The formulation of attractiveness equation is given in (2) 


B = Boe” (2) 


The distance between two fireflies is calculated using (3) 


k=n 
Si— sil = py (Sik — Sik)” 
= (3) 


Therefore, the final step in producing the new solution is done through (4) 





ri=| 


si= sit Boe” (s;— si) (4) 


From (1) and (2) the FA is controlled by the attractiveness, 6 and the fixed absorption coefficient, y. 
If the y of the current position for i-th firefly is more than j-th firefly, then it will be 6 = fo, Thus, it will move 
via attractiveness using (4). Otherwise, the firefly will move randomly. The general processes of FA are 
shown in Figure 2. 


2.2. Proposed Problem Formulation 

In this paper, the control variables in the FA optimization engine are the capacity of DG in the form 
of injected reactive power, Qpg. The reactive power is limited between Qpg min and Qpc max, i.e. the 
minimum and maximum capacity of the DG. The objective function is to minimize the real power loss in the 
system. The analysis will be conducted in IEEE 30-bus test system, hence there are thirty locations to be 
determined for DG placement. The overall optimization problem can be summarized as follows: 


Minimize { f(Qpg1, Qpa2,---Qpai,... Qoen) = Pross } (5) 
Subject to 

0 p.u. < FVSImay < 0.7 p.u. (6) 
0.85 p.u. < Vinin < 1.05 p.u. (7) 
Opa_min< Ong < QpGg_max (8) 
UPg = UPp + Prioss (9) 
2XOG = XOp + Qioss (10) 


Control variables in the FA are 


Si = [ Qne1, Qpc2,---Qnei,... Oven] (1) 
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Figure 2. Flowchart of FA 


Where, 

FVSInax : fast voltage stability index 

Vinin : Minimum voltage magnitude in per-unit 
Pross : Real power losses 

Qhoss : Reactive power losses 

Ong : DG’s injected reactive power 

xPG : Total real power generation in the system 
xPp : Total real power consumption 

XOG : Total reactive power injection in the system 
XOp : Total reactive power consumption 

Si : i-th individual in Firefly Algorithm 


Constraint (6) is significant for overall system stability assessment, while constraint (7) ensures that 
the voltage does not violate the system requirement. In this study, the lower limit for the voltage is set to be 
0.85 per-unit rather than 0.95 per-unit due to the increasing reactive loading that will be used in the 
experiment. Constraint (9) and (10) is to guarantee the power balance between the point of injection and 
extraction (i.e. generation and consumption) for both real and reactive power. 


2.3. Loss Sensitivity 

Prior to sizing the DG’s injected reactive power, the Loss Sensitivity (LS) technique will be applied 
for location identification in DG placement problem. The method is based on the change in losses with 
respect to generation increment [26], or in this case, the DG’s injected power. The terms (12) and (13) 
represent the loss sensitivity at k-th bus in terms of real and reactive power respectively. 


a _ pb 
L S, = Loss Pros 
a (12) 
a b 
15, = Yiu =o 
AQ, (13) 
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Where, 
LS; : Loss sensitivity at k-th bus in a power system 
b b 
Fross > Qioss : Real and reactive power losses before generation increment 
a a 
Fross » Qioss : Real and reactive power losses after generation increment 
AP, , AOx : Generation increment or DG’s injected power 


To obtain the losses before and after increment, at first, a load flow simulation is run for the first 
time to obtain the base case losses. Then, the generation increment is set at k-th selected bus and followed by 
another load flow simulation to obtain the new losses. The process is repeated until all buses has been tested 
and a rank of buses with corresponding LS indices is tabulated. 


3. RESULTS AND DISCUSSION 

Determining the suitable locations in IEEE 30-bus test system for DG placement was based on the 
Loss Sensitivity (LS) technique. In order to calculate the LS indices at all buses, the DG’s injected reactive 
power is set to a fixed value, in this paper, 10 MVar. Based from the analysis, there are three buses with 
significant LS indices, which are bus 30, 29 and 27 in descending order. Table 1 tabulates the results. 
Based on Table 1, there are two DG installations to be considered: double and triple installation. As has been 
found in Table 1, double installation involves placement of DG at bus 29 and 30, while triple installation will 
be at those two buses plus bus 27. The performance is measured in terms of minimum voltage magnitude, 
real power losses and overall system stability with respect to reactive loadings. The overall system stability is 
measured via Fast Voltage Stability Index (FVSI) [27] as the indicator. 

Table 2 tabulates the results when the reactive loadings were imposed on bus 30. All the 
improvement on voltage, losses and stability with respect to reactive loadings are illustrated in Figure 3 to 
Figure 5. The verification on the effectiveness of the proposed FA-based DG placement was conducted with 
bus 30 as the test bus for reactive loading. 


Table 1. LS Indices for Top 10 Buses in IEEE 30-Bus Test System — Injected Qpg = 10 MVar 








Rank Bus LS indices 
1 30 0.21674 
2 29 0.14861 
3 27 0.09791 
4 25, 0.07990 
5 26 0.06113 
6 24 0.04886 
7 23 0.02961 
8 22 0.02593 
9 21 0.02545 
10 28 0.02240 





Table 2. Effects of DG Installation on Voltage, Losses and Stability — Reactive Loading at bus 30 








oay Qpai Qpa2 Qpq3 Vinin P Loss F VSI, max 
Ooj.: sConciion (MVar)  (MVar) — (MVar) (p.u.) (MW) (p.u.) 

Without DG F : : 0.7424 20.7998 0.4918 

20 ‘Triple DG 11.6290 2.6661 41.3647 0.9908 18.2178 0.2555 
Double DG 18.7585 20.1940 : 0.9838 18.3193 0.2652 
Without DG E Ey - 0.7322 21.0002 0.7881 

30 ‘Triple DG 14.8401 12.4919 ‘12,0949 0.9402 18.2450 0.2184 
Double DG 41.0217 16.7647 = 0.9878 18.7495 0.3038 
Without DG : 2 . 0.7231 22.2271 0.7995 

40 Triple DG 18.8929 33.9285 ‘10.2528 0.9858 19.1583 0.3419 
Double DG 41.0217 16.7647 Z 0.9878 18.7495 0.3038 
Without DG “ : : 0.7100 23.2322 0.8332 

50 ‘Triple DG 47.6285 2.0685 13.4324 0.9681 17.8224 0.1940 
Double DG 75.7127 1.8545 0.9875 19.0330 0.4308 





Improvement on the minimum voltage magnitude seems to be consistent as shown in Figure 3. 
This is true for double DG installation where the voltage has been improved to a satisfactory level at all 
loadings. For instance, the minimum voltage has been improved from 0.7424 per-unit to 0.9838 per-unit at 20 
MVar loading, 0.7322 per-unit to 0.9878 per-unit at 30 MVar loading, 0.7231 per-unit to 0.9878 per-unit at 
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40 MVar loading, and 0.7100 per-unit to 0.9875 per-unit at 50 MVar loading. As has been expected, the real 
power loss is also been minimized through double DG installation as shown in Figure 4. At every reactive 
loading, it manages to keep the real power losses below the pre-optimization values, i.e. the values when 
there is no DG installed into the system. This is proven as follows: the double DG installation results in 12 
percent losses reduction at 20 MVar loading (i.e from 20.7998 MW to 18.3193 MW), 11 percent losses 
reduction at 30 MVar loading (i.e from 21.0002 MW to 18.7495 MW), 16 percent losses reduction at 40 
MVar loading (i.e from 22.2221 MW to 18.7495 MW), and 18 percent losses reduction at 50 MVar loading 
(i.e from 23.2322 MW to 19.0330 MW). 


1 24 
—6— Without DG 
0.95 23, |—*—Triple DG 
—®— Double DG 
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Figure 5. Effect of DG installation on stability-reactive loading bus 30 


Still, the reduction percentage signifies a good performance from the DG placement. The overall 
system stability, as shown in Figure 5, has been successfully enhanced through the reduction of FVSI index. 
The triple DG installation manages to keep the value below 0.4 per-unit regardless of the loading. 
For example, the FVSI,,,, has been reduced from 0.4918 per-unit to 0.2555 per-unit at 20 MVar loading, 
0.7881 per-unit to 0.2184 per-unit at 30 MVar loading, 0.7995 per-unit to 0.3419 per-unit at 40 MVar 
loading, and 0.8332 per-unit to 0.1940 per-unit at 50 MVar loading. From this analysis, it is not a mistake to 
say that the proposed DG placement using Firefly Algorithm (FA) is a good approach for any system 
operator (SO) of utilities to enhance distribution network performance, either in terms of stability, reliability 
and energy saving. This analysis justifies two findings: (1) the DG locations (i.e. bus 27, 29 and 30) are the 
most suitable locations for placement and (2) the DG’s MVar capacity (i.e. Qpg) as optimized by Firefly 
Algorithm (FA) is optimal. 
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4. CONCLUSION 

In conclusion, this paper has presented a method for optimal sizing of distributed generation (DG) 
using Firefly Algorithm (FA) and Loss Sensitivity (LS) technique for voltage stability enhancement. 
The optimization technique used in finding the optimal DG size is FA, while the LS technique was used as 
the indicator for the most suitable locations of DG placement. In the test system, two types of placement have 
been conducted namely double and triple DG installation. This study has revealed that FA-based DG 
placement, either double or triple installation, promises satisfactory results as the voltage magnitude, real 
power losses and system stability have been successfully improved. Hence, it can be said that the proposed 
method has solved the optimal DG allocation problem in distribution network. Therefore, it is aspired that the 
proposed FA-based DG placement will be extended to take into consideration some aspects ignored in this 
study, such as renewable energy sources (RES), investment cost, societal and environmental impact. 
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